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. Abstract: In this work we present a new approach on studying dynamical systems. 

Combining the two ways of expressing the uncertainty, using probabilistic theory and cred- 
ibility theory, we have research the generalized fractional hybrid equations. We have 
^ ^ introduced the concepts of generalized fractional Wiener process, generalized fractional 

Q ■ Liu process and the combination between those two, generalized fractional hybrid process. 

CN . Corresponding generalized fractional stochastic, respectively fuzzy, respectively hybrid dy- 

^ . namical systems were defined. We applied the theory for generalized fractional hybrid 

i Hamilton-Pontryagin (HP) equation, generalized fractional Hamiltonian equations. From 

O ' the general fractional hybrid Hamiltonian equations, fractional Langevin equations were 

O 



found and numerical simulations were done. 
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1 INTRODUCTION 

Fractional theory has applicability in many science fields. This approach presents frac- 
tional derivatives, fractional integrals, of any real or complex order. Fractional calculus 
is used when fractional integration is needed. It is used for studying simple dynamical 
systems, but it also describes complex physical systems. For example, applications of 
the fractional calculus can be found in chaotic dynamics, control theory, stochastic mod- 
eling, but also in finance, hydrology, biophysics, physics, astrophysics, cosmology and 
so on ([U], [in], [H], [12] )• But some other fields have just started to study problems 
from fractional point of view. In great fashion is the study of fractional problems of 
the calculus of variations and Euler-Lagrange type equations. There were found Euler- 
Lagrange equations with fractional derivatives, and then Klimek found Euler-Lagrange 
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equations, but with symmetric fractional derivatives [18j. Most famous fractional in- 
tegral are Riemann-Liouville, Caputo, Grunwald-Letnikov and most frequently used is 
Riemann-Liouville fractional integral. The study of Euler-Lagrange fractional equations 
was continued by Agrawal [2] that described these equations using the left, respectively 
right fractional derivatives in the Riemann-Liouville sense. This fractional calculus has 
some great problems, such as presence of non-local fractional differential operators, or the 
adjoint fractional operator that describes the dynamics is not the negative of itself, or 
mathematical calculus may be very hard because of the complicated Leibniz rule, or the 
absence of chain rule, and so on. After O.P. Agrawal's formulation [2] of Euler-Lagrange 
fractional equations, Baleanu and Avkar ^ used them in formulating problems with La- 
grangians linear in velocities. Standard multi-variable variational calculus has also some 
limitations. But in [31] C. Udri§te and D. Opri§ showed that these limitations can by 
broken using the multi-linear control theory. 

Another aspect that we use is the stochastic approach. Stochastic concepts were firstly 
introduced by J.M. Bismut, in his work from 1981, when stochastic Hamiltonian system 
was introduced. Since then, there has been a need in finding out tools and algorithms 
for the study of this kind of systems with uncertainty. Bismut 's work was continued 
by Lazaro-Cami and Ortega ([21], [22]), in the sense that his work was generalized to 
manifolds, stochastic Hamiltonian systems on manifolds extremize a stochastic action 
on the space of manifold valued semimartingales, the reduction of stochastic Hamiltonian 
system on cotangent bundle of a Lie group, a counter example for the converse of Bismut's 
original theorem. 

A new way for expressing the uncertainty is given in credibility theory. In this case we 
are not working on a probability space, like in the stochastic case, but on a credibility one. 
Credibility theory is based on five axioms from which the notion of credibility measure 
is defined, and it was introduced in order to measure a fuzzy event. This was first given 
by Li and Liu in their work p3]. This is a new theory that deals with fuzzy phenomena. 
Fuzzy random theory and random fuzzy theory can be seen as an extensions of credibility 
theory. A fuzzy random variable can be seen as a function from a probability space to the 
set of fuzzy variables, and a random fuzzy variable is a function from a credibility space 
to the set of random variables [23]. In our actual research, we will use fuzzy differential 
equations, that were firstly proposed by Liu [21]. This is a type of differential equation, 
driven by a Liu process, just like a stochastic process is described by a Brownian motion. 

In the case when fuzziness and randomness simultaneously appear in a system, we will 
talk about hybrid process. In this sense, we have the concept of fuzzy random variable 
was introduced by Kwakernaak ([IS], [20])- A random variable is a random variable 
that takes fuzzy variable values. More generally, hybrid variable was proposed by Liu 
|25j to describe the phenomena with fuzziness and randomness. Based on the hybrid 
process, we will work with differential equations characterized by Wiener-Liu process. 
This can computed using Ito-Liu formula [32] ■ In some situations, there exist many 
Brownian motion (Wiener process) and Liu process in a system, therefore, we can take 
into consideration also multi-dimensional Ito-Liu formula. 

In this paper, we restrict our attention to stochastic fractional Hamiltonian systems 
characterized by Wiener processes and assume that the space of admissible curves in 
configuration space is of class C^. Random effects appear in the balance of momentum 
equations, as white noise, that is why we may consider randomly perturbed mechanical 
systems. It should be mentioned that the ideas in this paper can be readily extended to 
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stochastic Hamiltonian systems [26j driven by more general semimartingales, but for the 
sake of clarity we restrict to Wiener processes. 

The paper is organized as follows. In Section 2 we present generalization of Riemann- 
Liouville fractional integral, Wiener process, and we have defined the generalized fractional 
stochastic equations. To get to a hybrid process, we have defined a generalized Liu pro- 
cess and the generalization of fractional fuzzy equation. The mixture between generalized 
Wiener process and generalized Liu process results as the generalization of hybrid frac- 
tional differential equations. In the third section we used the notion presented in Section 2 
for defining the generalized fractional hybrid HP equations. We have defined the general- 
ized fractional Riemann-Liouville, respective Ito, respective Liu integrals and in Theorem 
1 we gave generalized fractional hybrid HP equations. We have also defined generalized 
fractional hybrid Hamiltonian equations. The first order Euler scheme is presented and 
implemented for particular parameters. 



2 GENERALIZED FRACTIONAL HYBRID EQUA- 
TIONS 

Let / : M ^ M be an integrable function, a : M ^ M a function, p G M, p > 0. 
A Riemann-Liouville generalized fractional integral [12j is defined by 

ujtm = f / ,,j m - s)"(^-*)-ie-^(^-*)d., (1) 

ito r(a(s - t)) 

and 

T^{a{s-t))=T{a{z)%=,_t, (2) 
and T{a{z)) is Euler Gamma function given by 

T{a{z)) = {s- t)<'^-^e-^'-'^dt. (3) 



If a{z) = a = const, < a < 1, p = 0, from ([T]), results that 



1 '•^ 



^oItf{t) = Y^J^ f{s){t-sr-'ds. (4) 



Formula (jlj) is the fractional Riemann-Liouville integral, |15j . 

Generalized fractional Riemann-Liouville is a mixture between a fractal action used 
in physical theory and discount action with rate p, [12j. 

In the relations ([T]) and (Hj), s is called intrinsic time and t is called observed time, 
t ^ s. 

From ([T]), results that 



to' 

where 



It fit) = [ f{s)gt{s)ds, (5) 

J to 



I i{a[s - 1)) 
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Let {Q, 3^, Pr) be a probabilistic space characterized by the usual conditions, and 
(W(t))teR a 1-dimensional Wiener process. 

It is called a generalized fractional Wiener process, the process 

J-(f)= f g1{s)dW[s), t^s, (7) 

J to 

where is the function given in 

If x{t) = x{t,Lj) is a stochastic n-dimensional process and a : M x — M", b : 
M X M" ^ M", are deterministic functions, we will denote by 



tj^a{t,x{t)) = Jla{s,x{s))gf{s)ds, 
toJ?h{t,x{t)) = j!b{s,x{sM{s)dW{t 



(8) 



where a : M — > M, /5 : M ^ M, functions, the generalized fractional Riemann-Liouville 
integral, respectively generalized fractional Ito integral. 

We call generalized fractional stochastic differential equation, the functional Volterra 
type equation given by 

x(t) = x{tQ)+ j a{s,x{s))gt{s)ds+ [ b{s,x{s))g^{s)dW{s). (9) 

Jto Jto 

Using the notations given (IHl), it results that 

x{t) = x{to) I^ait, x{t)) +t, J^bit, x{t)). (10) 
The equation (IIUI) can be written formally in the following way 

dx = a{s, x{s))g"{s)ds + b{s, x{s))g^{s)dW{s). (11) 

Let us consider a(t,x(t)) = fi{t)x{t), b{t,x{t)) = a{t)x{t), where : M ^ M, cr : M — >• 
M, X : M — > R and a{z) = ai, j3{z) = ^-^y^- Then equation becomes 

1 /"* Li(s)x(s) , 1 /"* a(s)x(s) , , , 

x(t) = X to + TTT-^ / \: ds + — / 7 ^ ',,^ ds. 12 

^ ^ ^ °^ r(ai) ^ (t - s)i-"i r(i^) ^ (t - s)(i-"i)/2 ^ ^ 

The equation fll2p is called fractional differential equations that governs the stock model 
(Black-Scholes), ([H, [10], [S], [27], [29]). 

Let (9, y, Cr) be the credibility space with the usual conditions and {Lt)tm. an 1- 
dimensional Liu process [21]. 

We call a generalized Liu process, the following process 

i^"(0= f gns)dL{s), t^s. (13) 

If x{t) = x{t, 9) is an n-dimensional fuzzy process and Oi : M x R" ^ R", 6i : R x R" — >• 
R", are deterministic functions, we will denote by 



toHtai{t,x{t)) = Jl_^ai{s,x{s))g^{s)d{s), 
toLlbi{t,x{t)) = jl^bi{s,x{s))g^t{s)dL{s), 
4 
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where a:]R— s>R, >M are functions, and (7", g^ are given by (E]), the generahzed 

Riemann-Liouville integral, respectively the generalized Liu integral. 

We call generalized fractional fuzzy differential equation, the functional Volterra type 
equation given by 



x{t) = x{to) + Jl ai{s, x{s))g^{s)ds + x{s))g'^{s)dL{i 
= xito) +t, H^ai{t,x{t)) +t, L^h{t,x{t)). 



(15) 



Equation 0151) can be written formally as 

dx = ai(s, x{s))g"{s)ds + bi{s, x{s))gf{s)dL{s). (16) 

If ai{t,x{t)) = lJ,{t)x{t), bi{t,x{t)) = cr{t)x{t), fi : R ^ R, a : R ^ R, x : R ^ 
i, a{z) = 1, P{z) = l3i, from (JTOl) results that 

* ' a(s)x(s) 



x{t) = x{to) + ^(^)^(^)^^ + dL{s). (17) 

The equation (fT7|) is called fuzzy equation of a stock model [29] . 

Let (0,y, Cr) X {Q,3^,Pr) be the chance space [21], with the usual conditions, and 
iWt)tm an 1-dimensional Wiener process and {Lt)t^^ an 1-dimensional Liu process. Let 
x{t) = x{t, io, 9) an n-dimensional hybrid process and 02 : Mx M" ^ R^, 62 : Kx M" ^ M", 
C2 : ffi X M" — > M", deterministic functions. 

It is called a generalized fractional hybrid differential equation, the functional Volterra 
type equation given by 

x{t) = x{to) + j'l 02(3, x{s))g^{s)ds + 62(5, x{s))g^{s)dW{s) ^^^^ 
+ llc2{s,x{s))gt'{s)dL{s), 

where a : M [0, 1], /3 : M ^ M, 7 : M ^ M, are functions. 

With the notations given in ([8]) and (fT3|) . equation (fT8l) can be written as 

x{t) = x{to) +t, /°a2(t, x{t)) +t, Jt%2{t, x{t)) +t, K^C2{t, x{t)). (19) 

Formally, the equation ( fTSl) can be expressed as 

dx = 02(5, x{s))gf{s)ds + 62(5, x{s))gt {s)dW{s) + 02(5, x{s))g]{s)dL{s). (20) 



3 GENERALIZED FRACTIONAL HYBRID HP EQUA- 
TION 

Let Q be the paracompact configuration manifold and J^{R,Q) = R x TQ, T*Q 
the associated bundle of Q. Let {Q,7,P) be a probability space and (W{t),3^t)t£[a,b], 
where [a, b] C R, W(t) is a real-valued Wiener process and 3^t is the filtration generated 
by the Wiener process [5j. The HP principle unifies the Hamiltonian and Lagrangian 
description of a mechanical system. The classical HP integral action will be perturbed 
using deterministic function 7 : Q — M. 
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Let L : J^(]R, Q) — > IR be a function, called Lagrangian for the mechanical system 
and 7i, 72 : Q ^ M two functions of class C^. It is called generalized fractional action of £, 
with respect to the process (W(t))teM and (L(t))tg]R, the function A"" : QxVtx C{PQ) M 
defined by 

A^{t, q, v,p) = fjLis, q{s),v{s))+ < p{s), J - v{s) >)gfds + J^' ^Ms))9?dW{s) 
+ jll2{q{s))g?dL{s). 

(21) 

The first integral in (12T|) is called generalized Riemann-Liouville fractional integral, the 
second one is generalized fractional ltd integral and the third one is generalized fractional 
Liu integral. Moreover, 

C(PQ) = {{t,q,v,p) e C\[a,b], J\R,Q)), q G C\[a,b],R^), q{a) = qa, q{h) = g,}, 

[a, 6] CM, qaAbeW. 

We make the following notations q(t,6,uj) = q{t), v(t,6,uj) = v(t), p{t,9,uj) =p(t). 

Let c = {q,v,p) G C{[a,b],qa,qb) be curves on J^(M.,Q) between qa and g^, and B = 
{q,v,p,6q,6v,6p) G C°([a,6], J^(M,Q) x Q)) such that 5q{a) = 5q{b) = 0, and q, 6q 

are of class C^. 

Let {q,v,p){-,e) G C(J^(]R, Q)) be a family of curves on J^(M.,Q) such that they are 
differentiable with respect to e. The differential of the action A"' is defined by 



dA'^iSq, Sv, Sp) = ^^"(^, 0, qit, e),vit, e),p{t, e)) 



e=0 



where 

f^n(i\ — _ 



Sq{t) = §-q{t,e) , 6q{a) = 6q{b), 



e=0 



(22) 



6vit) = ^Xt,e) , 5p(t) = ^p(t, e) 
Using (j2Tl) . by direct calculus, we get the following theorem. 



e=0 



Theorem 1 Let £ : J^(R,Q) ^ be a Lagrangian function with respect to t, q 
and V and the first order derivatives are Lipschitz functions with respect to t, q, v. Let 
7ii72 : Q ^ M 6e functions of class C^, and with the first order derivatives Lipschitz 
functions. Then, the curve c = {q,v,p) G C{J^(R,Q) x M") satisfies the generalized 
fractional hybrid HP equations a.s. 

dq^ = v^ds, 

dp^ = {§-pMs,t))ds + ^dW{s) + ^dL{s), (23) 



Pi Qyi ) ^ 1, n, t ^ s, 



where 



,, , d(a(s — t)), , , a(s — t) — 1 1 dTi(a(s — t)) 

his, t) = l^ln\t -s\ + ^ ^ + p - —- ^. 

ds s — t Li[a[s — t)) ds 

□ 
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From ( 1231) we have: 
(i) If a{z) = 1, p = 0, then 

(ig* = v''ds, 



dp^ = §ds + ^dW{s) + '-^dL{s), (24) 



Pi Qyi ; ^ 1, ... , n, t ^ 

(ii) If a{z) = a = const, < a < 1, p = 0, then 

dp. = {§- P.^)ds + ^dW{s) + ^dL{s), (25) 

Pi Qj^i 1^ ^,...,Tl, t ^ S. 

For 72 = 0, we get the fractional stochastic HP equations fU]. 
If £ : M ^ M is hyperregular, that means 

d(^t[Sh) ^ 0' f™"^ (El results the 

following proposition. 

Proposition 2 (Generalized fractional hybrid Hamiltonian equations) 
The equations ( l23l) are equivalent with the equations 



d<t = ^ds, 

dp^ = i-§-p.h{s,t))ds + ^dWis) + '-^dL{s), 

where 



(26) 



H = piv' - L{t,q,v), 

da(s — t), , , a(s — t) — l 1 dTi(a(s — t)) 

his, t) = '-Int -s + ^ L + p - —- -T^^-^ ^ 

ds s — t Ti{a{s — t)) ds 

□ 

The equations fl2Ul) are called generalized fractional hybrid Langevin equations and can 
be written to describe the movement equations for relativistic particles with white noise 
and Liu process. 

Proposition Z If L = ]^gijV^v^ , where gij are the components of a metric on a manifold 
Q, then the equations ( l23l) take the form 

dq^ = v^ds, 

(27) 

dv' = -{Tij^v^v^ - h{s,ty)ds + g'^'^^dW{s) + g'^^^dL{s), i,j = 1, ...,n, 

where r*^^, are Cristofjel coefficients associated to the considered metric and h{s,t) is given 
above. 
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The equations ( l26l) become 

dPi = ii^P'p' - h{s,t)p,)ds + ^dW{s) + ^dL{s), I = 1, ...,n. 



Proposition 4 Let L : J^(M,M) ^ M 6e given by 

and 7i, 72 : M — > M. The equations f l28p are given by 
dq = pds, 



(28) 



□ 



(29) 



dp = - h{s, t)p)ds + '-^dW{s) + ™rfL(.). 

□ 

If V{q) = cos(g), 71(g) = ai sin(g) and 72(g) = |a2g^, the first order Euler scheme for 
the equations (!29l) is given by 

q{n + 1) = q{n) + Kp{n), 

(30) 

p{n + 1) = p{n) + K{sm{q{n)) — h{nK, t)) + ai cos{q{n))G{n) + a2q{n)L{n, Z2), 

where n = 0, iV — 1, K = ^, G{n) and L(n, z) are the simulations of Wiener and Liu 
processes and 

ainK — t) — I 1 

hinK, t) = a(nK, t)ln\t - nK\ + ^ — hp- ^ , , —VAainK - t)), 

riK — t Ti[a[nK — t)) 

with 

. t/tt(g - 1) ■ dVi{a{s-t)) 
a{s,t) = — , Ti{a{s,t)) = — , 

and 

G{n) = rando'm[normald[0,\/h]{l), 

2 "^^ 
^^^'^) = 1 , e.N/(/^^v/65.(n)) ' ^2(^) = 5^(6(«2g(fc),^(A:))). 

fc=0 

Using Maple 13, for the values of the parameters, a = 0.6, t = 0.8, ai = 0.1, 02 = 
0.3, 2:2 = 15., we get the following orbits. 



10 20 30 ■ 40 - 50 60 70 SO 90 100 t).022 0.024 0.026" - 0.02K- 0.030 



10 20 30 40 50 60 70 KO 90 100 



Fzgl: {n,q{n,Z2,uj)) Fig2: {n,p{n, Z2,uj)) FigS: {q{n, Z2,uj),p{n, Z2,uj)) 

8 



Chi§ Oana, Opri§ Dumitru 



If we are in the classical case, with a = 1, then the above graphics become 



10 20 30 40 50 60 70 KO 90 100 



10 20 30 40 50 60 70 HI) 90 100 



0.0205 0.0210 0.0215 0.0220 0.0225 



Fig4: {n,q{n,Z2,uj)) Fig5: {n,p{n, Z2,uj)) Fig6: {q{n, Z2,uj),p{n, Z2,u)) 

For a = 0.6, t = 0.8, ai = 0., a2 = 0.3, Z2 = 15., we get the orbits given in figures 7, 
and 9, and if a = 1, we get the figures 10, 11 and 12. 



10 20 30 40 50 60 70 fiO 90 100 

Fig7: ^2)) 



10 20 30 40 50 60 70 SO 90 100 

Fig8: {n,p{n,Z2)) 



0.03 0.04 0.05 0.06 0.07 O.OK 

Fig9: {q{n, Z2),pin, Z2)) 



10 20 30 40 50 60 70 fiO 90 100 



10 20 30 40 50 60 70 SO 90 100 



0.022 0.023 0.024 



FiglO: {n,q{n,Z2)) Figll: {n,p{n,Z2)) Figl2: {q{n, Z2),p{n, Z2)) 

The orbits for {n, q{n,uj)), {n,p{n,uj)), {q{n,uj),p{n,uj)), for the values of the param- 
eters a = 0.6, t = 0.8, ai = 0.1, a2 = 0, Z2 = 15., are represented in figures 13, 14 and 
15. 



10 20 30 40 50 60 70 fiO 90 100 

FiglS: (n, g(n, uj)) 



■ 10 20 30 ■ 40_ - 50 60_ 70 SO 90 100 

FigU: {n,p{n,uj)) 
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For a = 1, the figures 13, 14 and 15 become 



10 20 30 40 50 60 70 HO 90 100 

Figl6: {n,q{n,u>)) 



10 20 30 40 50 60 70 HO 90 100 

Fig 17: {n,p{n,u>)) 



0.0205 0.0210 0.0215 0.0220 0.0225 

Figl8: {q{n,uj),p{n,uj)) 



4 CONCLUSIONS 



In this paper we present generalization of fractional Riemann-Liouville integral, Wiener 
process, and we have defined the generalized fractional stochastic, Liu and hybrid equa- 
tions. The mixture between generalized fractional Wiener process and generalized frac- 
tional Liu process results as the generalization of fractional hybrid differential equations. 
We defined generalized fractional hybrid HP equations and generalized fractional hybrid 
Hamiltonian equations. The first order Euler scheme is presented and implemented for 
particular parameters. In the future work, we will consider other problems that deal with 
stochastic fractional HP principle. 
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